Nonequilibrium dynamics and thermodynamics of a degenerate Fermi gas across a 

Feshbach resonance 
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We consider a two-species degenerate Fermi gas coupled by a diatomic Feshbach resonance. We 
show that the resulting superfluid can exhibit a form of coherent BEC-to-BCS oscillations in re- 
sponse to a nonadiabatic change in the system's parameters, such as for example a sudden shift 
in the position of the Feshbach resonance. In the narrow resonance limit, the resulting soliton-like 
collisionless dynamics can be calculated analytically. In equilibrium the thermodynamics can be 
accurately computed across the full range of BCS-BEC crossover, with corrections controlled by the 
ratio of the resonance width to the Fermi energy. 

PACS numbers: 03.75.Kk, 03.75.Ss, 67.90. +z 



Recent experimental advances to control interactions 
in trapped degenerate gases by tuning (via magnetic 
field) through a Feshbach resonance (FR), led to ob- 
servation of long-sought after molecular Bose-Einstein 
condensation in bosonic and fermionic atomic gases 
0, 0] . Such tunability allows studies of these systems in 
previously unexplored regimes, as exemplified by recent 
JILA and MIT experiments 0,0, that appear to observe 
paired fermionic condensates along the crossover between 
the BCS regime of weakly-paired, strongly overlapping 
Cooper pairs, and the BEC regime of tightly bound, 
weakly-interacting diatomic molecules. An even more 
exciting possibility, unavailable in other related (e.g., 
superconductors and superfluids in condensed matter) 
systems, is the nonadiabatic switching of system's pa- 
rameters, thereby allowing access to highly coherent and 
nonequilibrium quantum states of matter. For bosonic 
85 Rb atoms, this was recently spectacularly realized in 
experiments by Donley et al []J. Using short magnetic 
field pulses that briefly bring the system close to a nearby 
FR, they observed coherent oscillations in the atomic 
condensate that can be interpreted as Rabi oscillations 
between atomic and molecular condensates Q. 

Stimulated by these experimental advances, in this pa- 
per we study a zero-temperature collisionless dynamics 
of a two-species degenerate atomic Fermi gas near a FR, 
that can be tuned through the Fermi sea. The goal is 
to understand the evolution of the system following a 
nonadiabatic change in an externally-controlled system's 
parameter, such as the detuning of the FR, loq, relative 
to the Fermi energy, ep. 

Our main result is the demonstration that such a sys- 
tem exhibits an integrable soliton-like solutions, corre- 
sponding to collective coherent oscillations of the Fermi 
gas between BCS- and BEC-like paired superfluid states. 
This dynamics is a paired-fermions analog of the atomic- 
to-molecular condensate Rabi oscillations observed by 
Donley et al in trapped bosonic gases. In the narrow 
Feshbach resonance limit we calculate the frequency and 
amplitude of these oscillations and find their analytic 



form. Our work directly builds on the recent discovery 
by Barankov, Levitov, and Spivak [f| of integrable dy- 
namics in the BCS model following a sudden change in 
the negative s-wave scattering length. 

Complementing our study of dynamics, we also analyze 
the thermodynamics along the full range of the BCS-BEC 
crossover |y, uj ■ We show that for a FR that is narrow 
compared to the Fermi energy, and a background scatter- 
ing length af, s that is short compared to the interatomic 
spacing, n^ 1 / 3 , low-temperature thermodynamics can be 
accurately computed analytically 

Our starting point is the Hamiltonian [t| 
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describing fermionic atoms, created by a£ CT with momen- 
tum p, "spin" a =1,1, and kinetic energy e p = p 2 /2m, 
that are coupled to diatomic molecular (resonant) states 
created by 6J. The position and the width (molecular life- 
time) of the FR are respectively controlled by the bare 
detuning energy e and coupling g, the former easily ex- 
perimentally tunable by a magnetic field. We have ne- 
glected nonresonant atom-atom and molecule-molecule 
interactions that we expect near a FR to be subdomi- 
nant to the resonant scattering. 

In what follows we focus on molecules with zero center 
of mass momentum, bo, neglecting molecules b q ^o ex- 
cited above the molecular condensate. For equilibrium 
phenomena, this approximation is justified at low tem- 
peratures and weak interactions (small g) , for which con- 
densate depletion is small. However, for nonequilibrium 
dynamics of interest to us, the validity of this approxi- 
mation is a more delicate issue, as one might generically 
expect such bosonic finite momentum excitations to be 
induced by a nonadiabatic shift in the FR. Nevertheless, 
physically we expect that for an initially homogeneous 
condensate and weak interactions the dynamics will 
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be dominated by 60 (t) ■ However, on sufficiently long time 
scales the b q ^o excitations and particle collisions should 
decohere and damp out the collective BEC-BCS oscilla- 
tions studied here, allowing a slow relaxation to a new 
equilibrium state for the shifted FR. Determining this re- 
laxational dynamics is beyond the scope of the present 
work. We thus replace b q in Eq. by &o<5 q ,o = &£q,o- 
Expected macroscopic occupation ((b'b) 3> 1) of the 
molecular level eo allows us to neglect quantum molec- 
ular fluctuations and replace operators b(t) , b' (t) by the 
corresponding c- numbers b(t), b*(t). 

We now look for time dependent fermionic opera- 
tors in terms of reference, time-independent fermions 
dyr, a~ P i, related to a p ^{t) through the Bogoliubov am- 
plitudes u p (t), v p (t) by a pT = u* d pT +v p a i _ pl , a_ pl = 

u*a^pi — v p a' p ^ with a constraint u*u p + v*v p = 1, 
that ensures fermionic anticommutation relations. The 
Heisenberg dynamics is then encoded in the equations of 
motion for u p (t), v p (t), and b(t) 

■f> 1 9b* .p. . gb f . 

id t Up = ~e p Up + — = v p , id t Vp = e p v p + —= u p , (2) 
v V VV 



idtb = eob 



fJ 



(3) 



The dynamical evolution equations J5J and J2J pre- 
serve pair correlations of fermions. Accordingly, we 
choose the initial state to be of the BCS type, |f2(Q)) = 

lip ( m p(0) + w p(0) a pT°-pi) 1^' ^ s a resmt °f tne time 
evolution the fermion wave function preserves the same 
form with time dependent factors u p (t) and v p (t). 

Following Ref. [5j we utilize Anderson's spin anal- 
ogy for the BCS problem and look for a vector- 
representation of these equations in terms of the variables 

Sp = Sp -+- iSp = 2Vpiip, Sp — VpVp u p^p: (4) 

that satisfy |S p | 2 = S*S P + S^S^, = 1. Equations © 
become 

d t S p = 2ie p S p - b*S p: d t S* = -2ie p S* + ^=bS p , 



d t s; = ^=(b*s;-bSp), 



(5) 



and have a form of Bloch equations for a spin 
precessing in an effective (p-dependent) field 
(— Re 6, Imfo, — 2e p ), whose azimuthal dynamics 
is in turn self-consistently determined by molecular 
(b(t)) evolution, Eq. J3J. It is now straightforward to 
check that these equations can be solved by an ansatz 
(for the BCS limit, in the absence of molecules suggested 
in Ref. @|) 



S, 



p(t) = e luJt (-(2e p - uj)C p Q(t) + %C p (l(tj\ , 



s;(t) = c p n 2 (t)~D p: 



gb(t) 
W 



(6) 



provided that the real function Sl(t) satisfies 

ft 2 + (n 2 - A 2 .) (O 2 - A 2 .) = 0, (7) 
and that C p and D p satisfy 



Dl-l 



" ' - A 2 _A^, 2^ = (2e p - uj) 2 + A 2 _ + A 2 + . (8) 



r<2 



Here A + , A_, and uj are parameters characterizing the 
periodic instanton-like solution Q(t) expressible in terms 
an elliptic integral. It follows from Eq. J7J that b(t)/\/V 
oscillates between the values A_ / g and A + / g (we choose 
A_ < A + ) with a period of oscillations given by 



2 / A 2 
A+ V A 2 



(9) 



where K(m) is the complete elliptic integral of the first 
kind. 

A direct substitution of this ansatz into Eq. J2J, shows 
that the solution, Eq. © is compatible with the evolution 
of b(t) provided that the following two conditions are 
satisfied: 
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(2e p - uj) C p , 



C n 



(2tt)3 p 



(10) 



(11) 



Eqs. (|10|l . Ullfl determine uj and A + in terms of the ex- 
perimentally controlled (initial conditions) parameters eo 
and A_ . As we will see below, the first of these equations 
is a nonequlibrium generalization of the BCS gap equa- 
tion. The second one simply reflects the conservation of 
the total particle number, ^ = 0, with 



N = 2b*b + V 



d 3 p 
(2Vp 



(s; + 1) 



(12) 



We note that Eqs. © only determine C p and D p up 
to a p-dependent sign, that one would expect to be fixed 
by the initial fermion momentum distribution, encoded 
in |f2(0)). With the exception of a filled Fermi sea ini- 
tial condition, O(0) = A_ = 0, the solution encoded in 
C p and -Dp, Eq. ©, does not correspond to initial con- 
ditions (itp(O) and i>p(0)) characteristic of a ground state 
of the Hamiltonian © for any value of bare detuning eo . 
Nevertheless we fix the sign of C p so as to most closely 
match n p = \{Sp + 1) to the initial fermion momentum 
distribution. For a large positive detuning this corre- 
sponds to a Fermi-Dirac step function with discontinuity 
at the chemical potential fi, that separates the hole-like 
and particle-like states. Combining this criterion with 
Eqs. ©, we find 



2 sign(e p - fj) 



[(2e P 



■uj 



A+ 2 ] -4A 2 _A 2 



= , (13) 
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where, as in equilibrium problem, fi is implicitly de- 
termined by the conserved total particle number N, 
Eq. CE2J, that reduces to N = V J As (1 - D p ). 

Equation involves a linearly-divergent inte- 

gral [llj , arising from an unphysical aspect of the model 
Eq. Q), that the atomic modes with an arbitrarily 
large energy interact with the molecular ones with equal 
strength g. In a more realistic model, the momentum de- 
pendence of the coupling g would cutoff this divergence. 
As usual, our ignorance of this high energy physics can 
be buried in a (uv cutoff-dependent) relation between 
the parameter eo appearing in the Hamiltonian and the 
position of the physical Feshbach resonance, luq. 

To see this, we calculate the two-atom scattering am- 
plitude within the model Hamiltonian Eq. Q . It is com- 
pletely determined by the self-energy of the molecules, 
given by 



E(£) 



dui d 3 p 



2tt (2tt) 3 (E - c 

g 3 p= 2 
i — m 2 v E — q" 

47T 



e p + i0+) (w - e p + i0+) 
(14) 



d 3 p 



rn 



(2tt) 3 p 2 



that diverges in exactly the same way as the integral in 
Eq. l|lUfl . Since T,(E) enters the retarded molecular prop- 
agator Gb{E) = (E — eo — S(_E)) _ in combination with 
eo, we can trade in eo for the physical ( "renormalized" ) 
FR detuning loq according to uq = eo — g 2 J jj^j ■ 
This leads to the two-atom scattering amplitude / = 
— }= 2 y= , that is consistent with the generic form 

based on unitarity |l2j and is identical to that of the 
Fano- Anderson model (equivalent to our model at the 
two-body scattering level) 0. We see that the scat- 
tering phase changes by 7r as the energy changes from 
below to above the physical ("renormalized") FR, ljo, 
with finite width for positive detuning, luq > 0, controlled 
by 7 = g 2 m 3 / 2 /47r. In contrast for negative detuning, 
luq < 0, the scattering amplitude has a real pole at neg- 
ative energy (a bound state in the o pen channel), that 
corresponds to a real molecular state |l4|. 

With above renormalization of detuning, the nonequi- 
librium gap equation l|10|) becomes 



UJq — UJ = 



d 3 p 
(2tt) 3 



2m 

(2e p -ui)C p - 



(15) 



where all the integrals are now convergent. This together 
with the atom conservation condition, (|ll|l . the total 
atom number Eq. Ijl2(l (that determines /i), the detun- 
ing wo, and the initial molecular density A_, allows us 
to determine the condensate frequency lu, and the molec- 
ular density maximum A + , that controls the period of 
oscillations, in accordance with Eq. ©. Although the 
complete solution requires a numerical evaluation of the 
integrals, here we will focus on two analytically tractable 
regimes: (i) A + — A_ <C A + that corresponds to small 



amplitude oscillations of the condensate about the BCS- 
BEC ground state, and (ii) A_ <C A + , that corresponds 
to the evolution of the filled Fermi sea, following a large 
downward shift in detuning. 

Let us first consider the regime of small oscillations, 
A + - A_ < A + about a BCS-BEC ground state. Limit- 
ing our analysis tow > 0, we find that Eq. I|ll|) constrains 
to to be close to 2/i. This together with the condition 
Eq. JIU) gives (at A_ = A + = A, and fi = lo/2) 



loq — 2fi = 



d 3 p 
J2vf 



1 



1 



( £p - M ) 2 + A2 € p 



.(16) 



This coincides with the BCS-BEC gap equation, that can 
be derived in the equilibrium treatment of this problem 
fl5j |. It relates the condensate density A/g to the Fes- 
hbach resonance luq, with the chemical potential /i de- 
termined by the total atom number equation, Eq. Ijl2|l . 
Simple analysis of these equations shows that for a large 
positive detuning, ojq ^> ep, molecules are strongly sup- 
pressed, leading the chemical potential to "stick" to ep, 
and to a conventional atomic BCS ground state, with 
A(w ) w 8e- 2 e F e-^- 2e "ys 2 H^) ^( e ) _ mM/^Tr 2 

is the atomic density of states). In this far off-resonance 
BCS regime the accuracy of the mean-field treatment of 
the equilibrium problem is controlled by the ratio of the 
width of the FR to Fermi energy, namely by the dimen- 
sionless parameter 7 2 /cf @- 

As the detuning ujq is lowered toward and below 6f, the 
chemical potential begins to track the detuning, ^t(wo) ~ 
wo/2 — C% 2 i/(wo/2)], with atoms from states between 
6f and /Lt(wo) converting into Bose-condensed molecules. 
The density of these tightly bound molecules, that co- 
exist with BCS's Fermi sea determine the gap, which 
displays a rounded mean- field behavior A cqu u(uJo) ~ 

J 2 ^r7(e^/ 2 — fi(uj ) 3 / 2 ). In the g — > limit, lu crossing 
of ep is a genuine quantum transition, with an upper- 
critical dimension of d uc = 2, and is therefore mean-field 
in 3d J6J. A finite atom- molecule coupling g rounds 
the transition into a smooth crossover near tp, that for 
small g (i.e., narrow FR, -f 2 <§; e^) is therefore also ac- 
curately described by the mean-field theory summarized 
by Eas. (|12llt)|) . Clearly, no additional anomalies appear 
when the 2-body FR (wo = 0) is crossed, since by that 
point nearly all atoms are bound up into well-ordered 
Bose-condensed molecules. One interesting observation 
is that in this fi < BEC regime the remaining di- 
lute fcrmionic atoms are paired but are nondegenerate, 
and therefore realize a "strongly-coupled" BCS state 01 • 
This picture of the ground state (that can be easily ex- 
tended to a finite temperature[l5j) is qualitatively con- 
sistent with recent observations by Regal, et al.Q, which 
find that molecules appear at about 0.5 Gauss above the 
experimentally determined 2-body FR (in our interpre- 
tation corresponding to ujq « 2e^)- 
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The nonequilibrium solution, Eq. then describes 
small oscillations about this equilibrium BCS-BEC state, 
with the period T = n/A given by Eq. |(?JJ. Because u> ps 
2/.1 it is easy to see that for a sufficiently small oscillation 
amplitude, the atomic momentum distribution function 
n p (t) — i(5p(£) + l) only changes near the Fermi surface, 
e p = li- 

The other interesting limit, A_ — > 0, A+ = A 
describes oscillations between a Fermi sea (a normal 
"metal", with fi = e F ) of atoms and bosonic molecules, 
following a nonadiabatic shift of detuning from oo down 
to ojq. In this regime, Eq. (|TT]l . (fT~5f) reduce to 



9 

UJq-UJ = — 



<r_ 

4 



d 3 p 
(2^)3 
d 3 p 



<=p 



sign(e p -e F ) 1 



C P 2) 4 

sign (ep - e F ) 



(2tt) 3 



A 2 
4 



(17) 



Simple integration of the atom-number conservation 
(second) equation above gives 

A w g 2 u(uj/2)tan~ 1 [(2e F -uj)/A}. (18) 

The behavior then depends qualitatively on whether the 
detuning ujq is shifted to the BCS (u)q/2 >> ep) or to the 
BCS-BEC (0 < uj a /2 < ep) regime, i.e., above or below 
the Fermi surface. In the former case Eq. I|18l) reduces 
to 2ep — uj ps A 2 /g 2 v(uj/2) <C A, with an exponentially 

suppressed BCS-like gap A oc e ^ v ^f) given by the gap 
(first) equation in Eq. (fTTj) . Hence in this regime, for 
narrow FR, oj ps 2ep and oscillations are confined to the 
vicinity of the Fermi surface. 

In contrast, for the detuning shift below the Fermi 
surface, < ujq/2 < ep, it is now Eq. ljT%|) that deter- 
mines the molecular density. It gives A ps g 2 v{uj /2) /2 = 
7 V / tJ x 0(1), while the gap equation gives ujq — uj = 
g 2 v{ep) x 0(1) -C ujo, enforcing the molecular energy 
u> to stick to the detuning ujq. We note, that, although 
A ss 7v / ^o is much larger (scales as g 2 rather than expo- 
nentially suppressed in 1/g 2 ) than for the detuning into 
the BCS regime, it is much smaller than the equilibrium 
A e gui;(wo); found as the solution to the equilibrium gap 
equation Eq. (|16fl in this regime. This suppression of the 
condensate oscillations is due to energy conservation be- 
tween the atoms in the Fermi sea and the molecules. For 
a narrow FR resonance, 7 v / wo < w» < ff it is only a 
small fraction of atoms in the Fermi sea that are in the 
immediate vicinity (set by the resonance width 7-y/wo) 
of luq, that can resonantly bind into molecules. The re- 
sulting resonant atom-molecule interconversion leads to a 
narrow oscillatory depletion of the Fermi sea, illustrated 
in Fig^ with a period of oscillations that diverges as 
T oc log ^-^^ , in the limit of the "normal" Fermi 

sea (A_ — > 0) as the initial condition. 

Analysis of the gap and atom-conservation equations 
shows that the amplitude of the atom-molecule oscilla- 




FIG. 1: Atomic momentum distribution function n p dis- 
playing oscillations between a Fermi-sea and Bose-condensed 
molecules. The atomic depletion in n p of width A ps 7 v / c^o 
appears at the molecular energy determined by the detuning 
emo- 



tions vanishes as the molecular energy uj ~ ojq approaches 
zero. This is again enforced by the energy conservation 
that in the absence of other degrees of freedom (e.g., 
molecules above a condensate, b q ^) forbids conversion 
of molecules at negative energy ujq into atoms at positive 
energy e p . We leave analysis of the dynamics that in- 
corporates these additional degrees of freedom to future 
research. 
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